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We study a totally asymmetric simple exclusion process (TASEP) with one defect site, hopping 
rate g < 1, near the system boundary. Regarding our system as a pair of uniform TASEP's coupled 
through the defect, we study various methods to match a finite TASEP and an infinite one across 
a common boundary. Several approximation schemes are investigated. Utilizing the finite segment 
mean-field (FSMF) method, we set up a framework for computing the steady state current J as a 
function of the entry rate a and q. For the case where the defect is located at the entry site, we 
obtain an analytical expression for J (a, q) which is in good agreement with Monte Carlo simulation 
results. When the defect is located deeper in the bulk, we refined the scheme of MacDonald, et.al. 
[Biopolymers, 6, 1 (1968)] and find reasonably good fits to the density profiles before the defect site. 
We discuss the strengths and limitations of each method, as well as possible avenues for further 
studies. 

PACS numbers: 05.70.Ln, 87.15.Aad, 05.40.-a 



I. INTRODUCTION 



Since its inception nearly four decades ago, the totally asymmetric simple exclusion process (TASEP) 
P, S H, 13, S Hi has become a paradigmatic model in non-equilibrium statistical mechanics. Not only is 
it one of the few mathematically tractable models in this field, it displays a rich variety of behaviors and 
provides insight to a range of complex physical systems, e.g., interface growth fC, biopolymerization 
[Til [l3 . ITsj and traffic [13, tlBl • In its simplest form, TASEP consists of particles hopping uni- 
directionally and stochastically on a one-dimensional (ID) lattice with complete exclusion (each site 
accommodating no more than a single particle). The original model [l| was defined on a ring (periodic 
ID lattice) and, despite having a trivial steady state distribution, displays complex dynamical phenomena. 
In a TASEP with open boundaries, there are even richer phenomena. Coupled to an infinite reservoir, 
particles enter/leave the lattice with rate a/ (3 (relative to the hopping rate within the lattice). The 
stationary state distribution was found analytically through a matrix ansatz ^ and displays three distinct 
phases along with continuous and discontinuous transitions (Q] . 

Independent of Spitzer [l| , a more general version of the open TASEP was proposed [9| to model the 
translation process in protein synthesis. In a living cell, the genetic code in the DNA is transcribed into 
messenger RNA's (mRNA's), which are then used to synthesize proteins (strings of amino acids) by a 
process which closely resembles a TASEP. However, to model this biological process properly, at least 
two major generalizations are required. Referring the reader to existing literature [1, [13, [HI, [3, [H, [3 
for the details, we only state these differences here: 

(i) Each particle "covers" t>l sites (typically 12 [ol. [TtI. [Ti] ) . i.e., exclusion occurs at a distance (.. 

(ii) The hopping rates are inhomogeneous, i.e., the rate for a particle at site i to hop (provided site i + i 
is empty) is 7^ and is expected to depend on the codon at i. 

Even the seemingly simple modification in (i) is so serious that an exact steady state distribution remains 
illusive. Only Monte Carlo simulations and mean field theories provide good estimates of certain steady 
state properties [13, [i3|- Though many properties are qualitatively similar to the I — \ case, such as 
displaying three phases (maximal current, MC, and low/high density, LD/HD) in the thermodynamic 
limit, there are important quantitative differences. For example, the phase boundaries in the a- (3 phase 
diagram shift to 



i.e., MC prevails for ck,/3 > x, LD for a < min(/3,x), HD for j3 < min(Q!,x). The average overall density 
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p and current J are also modified. In MC, we have p = \ — x and J = . In LD, p is now al/ (1 + at) , 
where 

£ = ^-1. (2) 
For HD, p actually remains the same: 1 — /3. In all cases, the current is given by 

J{p) = p{l-p)/{t-lp). (3) 

Beyond these simple quantities, the profiles are affected by ^ > 1 quite seriously [l^. ITgI. [20j . 

Clearly, generalization (ii) is much more intractable. A further complication is that the genetic code 
is not a "random" sequence. Thus, it is unclear if the notion of quenched random averages Elm,!!! 
- so successful in the studies of spin glasses [23| - is even meaningful here. Nevertheless, from the point 
of view of physics, it is reasonable to ask what the effects of inhomogeneities are on an open TASEP 
with extended particles. Along these lines, there have been several studies using different methods, on a 
variety of systems. Examples include a single "defect" in an otherwise uniform TASEP [ill, [H, E HI, 
two defects [ill HI], a cluster of defects [H,!!!], as well as a fully inhomogeneous set {7^} that is dictated 
by real genetic sequences [l3,[H,[l3|- In this context, we re-examine the open TASEP with a single defect 
here, i.e., 7^ = 9^ \]^i^k = 1- In particular, the steady state current is naturally suppressed if g < 1, 
but further, simulation studies found that this suppression is not as severe when the defect is located 
near the entrance (typically k < 0(10)) or the exit. This phenomenon was coined the "edge effect" [IB]- 

In this article, we focus on understanding this effect better. Since exact solutions are not available, we 
consider several levels of approximations, providing increasingly accurate predictions for the currents and 
density profiles. We also discuss the effects of introducing £ > 1 particles into the system. The remainder 
of this paper is organized as follows. In the next section, we provide some details of our model and a 
brief summary of previous results. As our approximations consist of neglecting certain correlations, we 
regard them as different levels of "mean field theories." Two new levels are presented in Section Hill In 
Section [IV[ we end with a summary and outlook for future research. 



II. THE MODEL AND SIMULATION DETAILS 



Our model consists of a ID lattice of N sites, with open boundaries. Each site, labeled by i — 1, 2, N, 
is either occupied or vacant, so that a configuration is specified by the familiar set of occupation numbers 
{rii} , Tii — 1,0. However, unlike the standard lattice gas model, we have particles of size i, in the sense 
that a single particle always occupies (or "covers") £ consecutive sites. Therefore, strong correlations in 
{rii} necessarily appear; not all 2-^^ possibilities of {n^} are allowed. A further complication is that, in the 
most-often used specification [l3, [l§|, a particle must lie fully on the lattice on the left (i.e., occupying 
i = 1, ■■■,£) but it can "dangle beyond" the right (i.e., only the particle's left most site must lie within 
the lattice). As a result, the total number of holes on the lattice can vary even for a given, fixed number 
of particles. This complication can be ameliorated, however, if we choose the lattice to have N + £ sites 
and symmetrize the rules for entrance and exit. To conform with the notation of previous studies, we 
will avoid this route here. 

An alternative specification is to locate each particle by one of the £ sites, e.g., its left most site. With 
protein synthesis in mind, we follow [l3| and refer to this special site (on the particle) as the "reader." 
The motivation comes from the ribosome "reading" the next codon (and waiting for the arrival of the 
associated transfer RNA) before it can move onto the next codon. With this convention, we define r.i — 1 
if site i is occupied by a reader and = otherwise. Clearly, {r^} labels a configuration and, like {rii}, 
there are strong correlations. Choosing to locate the reader at the left end of a particle [ij, {rii} can 
be generated from {ri} by rii+j = for j € [0, £ — 1]. We will also use the reader position to locate the 
particle, so that ri — 1 will be used interchangeably with "A particle is located at site i." Finally, we 
define all sites beyond N to be free, so that a particle at the last £ sites is not hindered sterically by any 
others. 

Turning to the dynamic rules, it is easiest to state them in terms of r^. The motion of an interior 
particle is obvious; only the entry/exit rules need clarification. Coined "complete entry, incremental 
exit" in Ref. [11], these are: 

• ri 1 with rate a, provided r^ ~ 0, for fc e [1,^]; 

• n — \ ^ 7'i4-i = 1 with rate 7^, provided ri^g — 0, for i e [l, — £ — l] ; 
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FIG. 1: Sketch of a TASEP for particle size £ = 6 with a single slow site at position k, with rate q. 

• Ti — 1 ^ 7-j+i = 1 with rate ji, for i e [N ~ £, N ~ l\; 

• rjv = l— »fAr = with rate (3. 

In our simulations, we establish an array of N entries to represent the lattice sites, as well as an extra 
one {i = 0) for the reservoir. We use a random sequential updating scheme and keep track of the locations 
of readers. In one Monte Carlo step (MCS), we make M + 1 attempts to update, where M is the total 
number of particles on the lattice. As the 1 accounts for a particle in the reservoir to be chosen, there 
is an even chance for each particle to be updated once, as well as introducing a new particle into the 
system. A sketch of this process is shown in Fig. [TJ The lattice is initially empty and we discard the first 
2 X 10® MCS to ensure that the system has reached the steady state. A further 2 x 10® MCS are used 
for collecting measurements, each separated by 100 MCS in order to avoid temporal correlations. Unless 
otherwise noted, averaging over the 2 x 10^ measurements provides good statistics. Such steady state 
averages will be denoted by (...). We studied different system sizes between 200 and 1000, with most data 
taken from N = 1000. 

To characterize the state of the system, we monitor several observables. The most obvious is 

Pi ^ in) , (4) 

a quantity we will refer to as the reader density. Of course, pi is just the average number of particles in 
the system (i.e., ribosomes on the mRNA). Thus, the overall particle density J2i Pi upper bound 

of Another interesting variable, pi = {rii), referred to as the "coverage density", is the probability 
that site i is covered by a particle (regardless of the location of the reader). Of course, it carries the same 
information as p\, since the two are related through 

J P^-l^k=oP^-k (5) 
I Pi ^ Px - Px-1 + PUi 

(with the understanding pj = for i < 0). The overall coverage density, p = jj X^i Pi' reach unity 
and provides a good indication of how packed the system is. From pi, we can also access the profile for 
the vacancies (holes): 

P- = 1 - P. ■ (6) 

When £ = 1, the two density profiles are of course identical. As soon as ^ > 1, serious correlations 
appear [H, [H, [H] . 

A quantity of great importance to a biological system is the steady state level of a given protein. If 
we assume that the degradation rates are (approximately) constant under certain growth condition, then 
these levels are directly related to the protein production rates. In our model, such a rate is just the 
average particle current J, defined as the average number of particles exiting the system per unit time. 
At steady state, it is also the current measured across any section of the lattice. For simplicity and to 
ensure the best statistics, we count the total number of particles which enter the lattice over the entire 
measurement period. 

For our investigations here, we focus on one simple type of inhomogeneity: a single "slow" site (Fig. [T]) 
in an otherwise homogeneous lattice, i.e., a bottleneck along a smooth road. Locating the defect at site 
k, we have 



Ji^k = 1 and jk = q 



(7) 
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with q < 1 . A common approach to this type of problems is to study the lattice as two sublattices (left, 
sites 1 to k, denoted by L, and right, the rest, denoted by R) connected by q and having the same through 
current Jl = Jr. We are especially interested in the dependence of the current, denoted by J{q,k), on 
the parameters q and k. 

Previous studies located the defect far from the system boundaries, e.g., k « N/2 [H, There, 
it is sufhcient to regard both sublattices as infinite and to exploit the results of the single TASEP while 
matching L and R appropriately. Of course, the matching condition is not exactly known and the previous 
studies propose different approximation schemes. These approaches lead to tolerably good predictions for 
the average densities and currents. Here, we provide examples for the a = f3 = 1 case. In the most naive 
scheme (referred to as the "naive mean- field," NMF, approximation in the exact expression for the 

current, q {r^ (1 — nfc+^)), is replaced by q (rfc) (1 — nj^j^n). Despite the severity of this approximation, the 
result for the current 

r 9/[(l + g)(l + g^)] for9<lM /o^ 
'^^^^ = \ ^2 forg>l/V£ 

captures all the features of the system qualitatively. An alternative approach fl3| takes into account some 
of the correlations in (rknk+f ), leading to 



JsKL = (l - Q - - 2Qj / {Qe} (9) 

where Q = 2qi (l + qtj j (l + + 2q(^ . The agreement with data is much better than NMF, as shown 
explicitly in [20|. 

However, it is clear that neither scheme can provide any information on how the system is affected by 
the location of the defect, k. Now, simulations with I up to 12 showed a non-negligible increase |16i . [20l . [25| 
in J as the defect approaches the system boundaries, a phenomenon coined the "edge effect" in In 
the next section, we consider two more refined approaches. One is based on a mean-field theory proposed 
by MacDonald, Gibbs and Pipkin (MGP) 0, which we generalized to an inhomogeneous TASEP. The 
other, proposed by Chou consists of a cluster approximation (FSMF, finite-segment mean-field jllj'). 
which accounts for the "interaction" between the entry rate a and the defect q. Though both provide 
much improvement over the expression in Eqn.® above, their limitations will be discussed. 



III. BEYOND SIMPLE MEAN FIELD THEORY 



To appreciate the various levels of approximations, we begin with the exact expressions for the current. 
From the dynamic rules, J is given generally by 

J = a (1 - ni) = ap\ (10) 

= 7, (r, (1 - ; ^e [l,iV-^] (11) 

= ; ze [7V-f,7V-l] (12) 

= /3(rAr) . (13) 

If we define = /3, then the last two equations are just J — ^ip\ for the last I sites. 
In the study here, we have ^i^k = 1 and small k. So, Eqn. (fTTI) becomes 

J = (n(l-n,+,)) ; ie[l,N-e],i^k (14) 

= q (rk (1 - rik+e)) ■ (15) 



A. An approach using the MGP recursion relation 

In the naive mean field approach, {riUj) is replaced by (r^) {rij). For ^ = 1, this approximation turns 
out to be quite good. Thus, the profile is well described by the solution to the one-term recursion relation, 
namely Jnmf = Pi (1 ~ Pi+i)- Unfortunately, this naive approach fails for £ > 1. The difficulty is due 
in part to (r^) 7^ (n^) and in part to the severe exclusion at ^ > 1. MGP took into account some of this 
exclusion 3 and proposed a much better approximation. The key lies in replacing p ( /i^ | or /i^), the 
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conditional probability of finding a hole at site k given that this site is occupied by either a reader or a 
hole, according to the fraction: 

p{hk \rk or hk) = ,. ' ^ ■ (16) 
Pi+Pi 

Notice that, for £ = 1, the denominator is simply unity and an equality holds. Now, (ri(l — rii^g)) is 
given by the product of pi and p {hi+e\ri), the conditional probability of having a hole at site i + £ given 
there is a reader at site i. But, if a reader exists at site i, then site i + £ must be occupied by either a 
hole or a reader, so that p {hi+i\ri) = p (hi+i \ri+i or hi+i). With the key approximation above, MGP's 
scheme can be summarized by 

(r,;(l - n^+i)) = p\p{h,+t\ri) /'^''^'^y, ■ (17) 

Pi+e + Pi+e 

Generalizing to the inhomogeneous case, an approximate Eqn. (jlip now reads: 

Jmgp =7. /'"j^ ■ (18) 

Pi+e + Pi+e 

To proceed, we follow MGP and regard this as an £-term recursion relation. Starting from the last site, 
Eqns. p2ll3p allow us to write down the first £ terms: 

= Jmgp/7»; « e [N-£,N] . (19) 
With these, we can start a "backwards" recursion (BR) relation: 



'^MGP L ^ Pi+e 



Jmgp 


1 J2k=l Pi+k 


It 


1 X]fe=i Pi+k 





(20) 



and obtain the rest of the densities (z = 1, N — £). Finally, to fix the unknown JmgPi we impose 
Eqn. HOl): 

JMGP-ajl-^p'fej (21) 

Though Jmgp is "just the solution to a polynomial equation," its exact value is quite intractable, since the 
order of the polynomial approaches 2^ for large £ (such as £ ~ 12). Unfortunately, numerical techniques 
are also of limited value, due to the extreme sensitivity of the BR to small inaccuracies. As a result, 
given the computational power of four decades ago, MGP were able to exploit this approach only for a 
limited range of £ and N. Our interest here is the edge effect associated with just one slow site near the 
entrance. So, we would be considering a short sublattice (L, length k < 30) coupled to a longer one (R). 
Thus, we are in an ideal position to exploit MGP's approach for the L sublattice, while matching it to 
the results of an infinite TASEP for the R sublattice. 

Noting that the particle densities are uniform after the slow site (i.e., pl^f^ = P]i = PrI^ ): we approx- 
imate the R sublattice as an infinite system, so that Eqn. ^ for a homogeneous TASEP applies: 

J _ Pi?.(l - PR) ^ Pfl(l - lp\) .22^ 

£-£pR l-ipR 

Of course, neither nor p^?, is known, and both must be determined through matching conditions to 
the L sublattice and q. For the latter, we first consider expression (fT8)) for site k: 



Jmgp — IPk ^ , h ^ IPk ' 



^P'r + p'r l-^P\ 

where we have inserted the uniform density noted at sites beyond k. Of course this must be J/j, which 
provides us with the reader density at the slow site: 



Pk = p'r/i- 



(23) 
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Continuing with the BR, we have 



Jmgp { 1 



Pk-i+e 



I, Z^n=l Pk-i+n 



Pr^-^Pr) i-ELipI-,+„ 



(24) 



(25) 



until we have all the reader densities, p\,...,p'f., as functions of p'^. Finally, we impose Eqn. (PT|) for our 
case 



1 _ PRi^-iPR) 

a 1 - £p'j^ 



n=l 



(26) 

12, A; = 26 and g = 0.2, and 
(27) 



which fixes the unknown p^?, ^ind so, all quantities of interest. 

As an illustration, we carry out this program for the specific case of £ 
find [11 

p\ ^ 0.0506; Jmgp = 0.0448 . 

These values agree reasonably well with the simulation results of (0.0544, 0.0472), respectively. The 
density profile pi associated with this result, shown in Fig. [21 is labled RRl and by the circles (red 
online). Apart from the first few sites, the fit is quite respectable. However, the fit improves considerably 
if we arbitrarily relax the constraint (j26p . In the same figure, we also display such an alternative (x's, 
blue online, labled RR2), obtained by choosing 



PR 



0.0560 



Jmgp = 0.0478 



(28) 



We found the coefficient of this fit improves from 0.97 to 0.99. Moreover, it provides much better 
agreements with the i?-sublattice density as well as the overall current. The price, however, is a rather 
poor a, with the right hand side of Eqn. ([26|l missing unity by 33%! It is unclear why the BR displays this 
peculiarity, although we should perhaps not expect much better agreements, given that some correlations 
are ignored in this approach. 
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FIG. 2: Density profile obtained through a BR relation, q = 0.2, k = 26 and £ — 12. The parameters used for 
fitting are J and pr. 

To summarize, we present the results of the two BR's in Table ID If we impose the constraint ([26]) 
seriously, we see that the overall fits are tolerable. On the other hand, if we relax this constraint, then 
there is substantial improvement on all quantities except "a" . Clearly, this method, however unsystematic, 
manages to capture much of the details of the edge effect. Unfortunately, the BR relation fails to produce 
the long tails in the reader profiles (e.g.. Fig. 5 of [l^l)- Indeed, Eqn. fTS)) becomes very unstable after 
about 40 steps, setting a limit on the maximum k to which it can be applied. We believe there are 
inherent difficulties with this approach (beyond that of machine accuracy), but that is outside the scope 
of this paper. 
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TABLE I: Summary of fit parameters and a for the two recursion relation schemes. See text for the details. The 
percentage deviation from the simulation is included in the parenthesis. The last row shows the value of B? from 
the fit to the profiles in Fig. 2. 





simulation 


RRl (%) RR2 (%) 


PR 


0.0550 


0.0506 (8.00) 0.0560 (-1.82) 


J 


0.0472 


0.0448 (5.08) 0.0478 (-1.23) 


a 


1.0 


0.95175 (4.82) 1.5045 (50.45) 
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FIG. 3: Sketch of one slow site g at fc = 1. FSMF matches a two-site TASEP with the rest of the system. 



B. Finite-segment mean-field (FSMF) theory 

Given that we are interested in the edge effect, we can improve on the above method by accounting for 
the physics of the small L sublattice exactly. This approach follows the work of Chou and Lakatos 
in which the "finite-segment mean-field theory" was developed to understand quantitatively the effects 
of clustered defects. Here, we generalize this method to particles of size t > 1 and solve the full master 
equation explicitly for the L sublattice (for small fc). The key idea is to find the exact expression for the 
current for this small finite segment and then match it to the result of an infinite system (i.e., the R 
sublattice). The approximations appear only in the matching conditions and finite size, O (l/N), effects 
associated with the latter. We further consider the interplay between the defect rate q and the on-rate a. 
Although our results are based on this simplified model, understanding such interactions elucidates the 
effects of having a slow codon or a cluster of slow codons near the initiation codon of an mRNA, which 
is frequently observed in living organisms [13, [111 . 

For the L sublattice, the maximum dimension of the transition matrix is 2^'^+^' , which requires enormous 
amounts of computing time even for fc = 10 and £ < 12. Fortunately, for k < £, there can only be 
i(fc^ + 3k + 2£ + 2) allowed states in the k + £ segment. In particular, for fc = 1, the dimension of the 
transition matrix is £ + 3 which is sufficiently small for us to understand the "edge effect." We illustrate 
this method through a detailed account for the simplest case (fc = 1), showing the results for both £ = 1 
and £ > 1 (as there are subtle differences for the latter). 

For the case of fc = ^ = 1, we consider two sites, i.e., the entire L sublattice and one site for the R 
sublattice. Thus, we have only two rates: the entrance rate a and the rate for the slow site q. There are 
four possible configurations (ni — 0,1; n2 = 0,1), labeled by xo,i,2,3- For convenience, we list these in 
their binary sequence, namely xq corresponds to state (0,0), xi to (0,1), etc. The master equation for 
the evolution of Pi, the probability to find the system in Xi at time t is 

dtP,{t)^Y.^wiP,{t)^w)m)] (29) 
j 

where Wj is the rate of Xi to Xj . Referring to Fig. [3] and writing the right hand side of the above as a 
matrix W operating on a vector P(t), we write explicitly 



W = 



/ 


—a 


Pl 





\ 







-a - f3L 


q 







a 





-q 


Pl 


V 





a 





-h) 
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Here, Pl is an effective exit rate, which is to be fixed by matching. The stationary state distribution P* 
is easily found: 



P* = z- 



( PL/a \ 
1 

{a + f3L)lq 
V a/PL J 



where Z = 1 + /J^/a + 0//?^ + (a + is the normahzation factor. The steady state current Jl foUows 
readily: 

Jl = a (1 - ni) = q (m (1 - ^2)) - Pl (n2) (30) 

= (a + /3i)/Z (31) 

Meanwhile, on the R sublattice, we have 

Jr. = {n2 (1 - "3)) = ("3 (1 - ^^4)) = ••■ 

which becomes Jr = pr{1 — pn) for an infinite system in the low density phase. Matching this to the 
last equation in (|30p . we arrive at 



and so, 



/3l = 1 - Pfl 



Jl = Jr^ Pl (1 - Pl) 



Setting this equal to expression (|3ip , we find an equation for [3l ■ The final answer for the current in this 
approximation scheme is 

(X 

^FSMF(a,g;£= 1) = ——-{{l + a){q-a)+'R} (32) 
2(g + a) 

where 7?. = ^J (\ -\- a){q + a)[a'^ + a + q — 'iaq). 

We caution that this formula should be used with some care. Though not very transparent, this function 
monotonically increases with both q and a when both are small. However, beyond a line in the a-q unit 
square, it decreases back to zero. The maximum J on this line is precisely 0.25, at which point the system 
enters the MC phase. Substituting 0.25 into the left of Eqn. (I32p . we find the phase boundary: 

a{2a + 1) 

^-(")^ 4a^ + 2a-l ^^'^ 

beyond which [q > qc) the MC state prevails. 

To further appreciate the quality of this theory, we present the comparison between its predictions and 
simulation data in Fig. |4l Specifically, we show three typical scans through the a-q plane, sketched in the 
inset of Fig.[31 The upper-right corner represents the MC phase (where J — 0.25), with the solid (black 
online) curve being qc (a) in Eqn. ([33]). The color-coding in the inset matches the data in the main plot 
of Fig. m Here, the statistical error associated with the simulations is estimated to be around 0.01%. 
We observe excellent agreement (within 1%) between the data and the theory results. Needless to say, 
we can extend this approach, with some labor, to fc > 1. Since we doubt that the agreement would be 
significantly different, we believe there is no need to pursue this investigation further. Instead, we turn 
next to the more interesting systems with extended objects. 

Generalizing to the case with £ > 1, we need to account for two novel aspects. One important difference 
is the current matching condition: Instead of — (3^ — — Pr), we use the MGP result for ^ > 1 case 
and 

I-PR ^^(Wg)_ 

The other new item is that, since we expect some "period-^" structure, we could consider all k up to 
£ and still restrict ourselves to having only one particle in the L sublattice. Following the spirit of our 
analysis above, we also wish to account for the effects (on the L sublattice) due to a particle which 
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FIG. 4: Comparison between the simulations(solid symbols) and the FSMF approximations from Eg . 1321 (connected 
open symbols). Top (black online; □): q = 1; Middle (red online; •, Q)- = g; and Bottom (blue online; 
T, o): g = 0.2. In all cases, £ = 1 and = 1000. The inset follows the same color coding scheme. See text for 
details. 



just moved into the i? sublattice. With at most two particles in our finite segment, the exclusion "at a 
distance" means that we need to study a system with k 1 sites. Fortunately, the configuration space 
(for, say, the reader occupations {r^}) is still manageably small. Thus, there is just one 0-particle state, 
(k + €) 1-particle states, and fc (fc + 1) /2 2-particle states. We demonstrate the case for fc = 1 here. 
Let us label the 0-particle state by xq, the state with a reader in site i by xi^^-i {i = 1, + 1), and 
finally, the 2-particle state by xe-^.2■ To be pedantic, we show the explicit set of occupation numbers {rj} 
corresponding to these x's: 

a;o^(00 Oo) 

xi ^ (^0 Ol) 

a;2^(00 lo) 



(35) 



Xl 4^ 
Xi+i ^ 

xe+2 ^ 



1 

1 
1 



o) 
o) 
oil 



The advantage of this slightly peculiar labeling is that it reduces to the £ = 1 case easily. 

To find the current J{a,q,£) for £ > 1, we need to compute the new transition matrix 
configurations clear in our minds, we simply write: 



With the 



W: 







/3l 
-a - Pl 






/3l 
-Pl 




\ 



•• 

Pl 

-/?L Pl 



a 

V 





a 







(36) 
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Similar to the £ = 1 case, we find 



1 

P* = Z' 

[a + 
V a/h J 

to be the same, except for £ more entries of (a + in the middle. Thus, 

Z=l + pL/a + a/pL + + pL)/q + i{a + Pl)/I3l 

Meanwhile, we still have Jj^ — [a + (3l)/Z. Finally, matching Jl with — Pr{^ — Pb.)/ ~ ^Pi?) and 
using Eqn. we arrive at the solution for general £ (and k = 1 ): 



a 



q-a + 7^(l + a)-^ +£ q{l - a) + Uq {q + a)' 



2(9 + a) 



(1 + + ^ [l + £aq{q + aY 



(37) 



We have written Jfsmf in a form that clearly reduces to Eqn. ([5^ ior £ = 1. Given the shifted phase 
boundaries for £ > 1, the system enters the MC phase when: 



qc = 



ax{a + 1 - x) 



a{a + 1 - x) - x(l - x) 



(38) 



In the next section, we will summarize our findings along with some comparisons to Monte Carlo simu- 
lation data. 



IV. SUMMARY AND OUTLOOK 



We investigated how a single defect site near the lattice boundary (small k) influences the steady state 
properties of the system for both point particles and those of length £ > 1. The simplest "mean- field" 
approaches - NMF and SKL - are unsuitable, since both rely on matching two infinite TASEP's across 
a defect and cannot address the issue of fc-dependence. Instead, we considered two more sophisticated 
levels of mean-field methods, with complementary strengths and weaknesses. One method, first used by 
MacDonald, et. al. [9|] (MGP), is based on a recursion relation for the density profile. The advantage of 
MGP is that, up to moderate k values, its predictions for both the profile and the current are reasonably 
good. The weakness is that we can access these predictions only numerically so that the dependence on 
the control parameters a, q, k, and £ remains obscure. It is also unclear how to systematically improve on 
this approach. The other method, based on an exact account of the physics of the first k + £ sites, is a 
generalization of the finite segment mean-field (FSMF) theory of Chou and Lakatos [ll| . The strengths of 
this method are many. Based on the steady state solution to the full master equation, it can be improved 
systematically. Its predictions agree with simulations exceedingly well and provide analytic expressions 
so that the dependence on {a, q, £) can be appreciated. Both methods are obviously severely restricted to 
a relatively small range of fc's. In MGP, the limitation arises from the extreme sensitivity of the recursion 
relation and, in FSMF, an exponential increase (in the worst scenario) in the size of the transition matrix. 
Moreover, some of the long tails in the "edge effect" extend up to fc ~ 50 '^^^^ beyond the present 

reach of either approach. Hopefully, more efficient approaches will be developed in the future. 

To summarize, we find three successively better methods to describe a TASEP with a defect near the 
entrance, for t > 1. To illustrate, we present the results of all three mean-field approaches, along with 
simulation data, in Tables HIl and Hill Both concern the case with a — (3 = 1 and fc = 1; the difference 
being £ = 1, 2 in the two Tables. We see that, by accounting for exclusion at a distance, MGP (with the 
fit through RRl) clearly succeeds better than NMF when £ > 1. Meanwhile, it is hardly surprising that 
an exact treatment of the finite segment before the defect is superior to both. 

Though not displayed explicitly, similar improvements are found to hold for £ up to 12. Finally, if 
we use Jfsmf for k = 1 and Jnmf for fc ^ 1 (i.e., deep in the bulk), we arrive at a prediction for 
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TABLE II: Different mean-field approximations for the current J{q) for the case a — f} — 1, k — 1 and £ — 1. 
Simulation results are based on a lattice with L = 1000. 



9 


J'nmp 


Jmgp 


JpSMF 


simulation 


0.1 


0.0826 


0.0833 


0.0863 


0.0864 


0.2 


0.1389 


0.1421 


0.1498 


0.1490 


0.3 


0.1775 


0.1848 


0.1954 


0.1941 


0.4 


0.2041 


0.2153 


0.2261 


0.2248 


0.5 


0.2222 


0.2361 


0.2440 


0.2432 


0.6 


0.2344 


0.2476 


0.2500 


0.2502 



TABLE III: Different mean-field approximations for the current J{q) for the case a = f3 = 1, k — 1, and £ = 2. 
Simulation results are based on a lattice with L — 1000. 



q 


Jnmp 


Jmgp 


JpSMF 


simulation 


0.1 


0.0758 


0.0763 


0.0788 


0.0771 


0.2 


0.1190 


0.1213 


0.1266 


0.1231 


0.3 


0.1442 


0.1484 


0.1543 


0.1509 


0.4 


0.1587 


0.1639 


0.1680 


0.1665 


0.5 


0.1667 


0.1708 


0.1715 


0.1717 



Ai (q) = Jq(k = 1)/ Jq{k — !■ oo), defined in (20| . The remarkable non-monotonic behavior in Ai (shown 
in Fig. 7 of [2^) is well captured by the combination of these two mean-field approaches. 

Beyond our investigations here, there is ample room for future research. In an open TASEP, there 
are two "edges" and so, two possible "edge effects." We reported findings for only one. When the slow 
site is near the exit {k ~ L), the current is also observed to increase [1^. However, due to lack of 
particle- hole symmetry for £ > 1, this increase is not the same as the case for small k. Further, there 
are serious complications associated with the profile, especially for small q (e.g.. Fig. 5 in [l^l)- Thus, 
it would be desirable to find better methods to understand these peculiarities quantitatively. Similarly, 
we should explore the "edge effect" for q > 1. When a "fast site" is deep in the bulk, it has little effect 
on the current. However, its effects if located near the edges, especially near the exit end, are yet to 
be discovered. Beyond one defect, there are obvious questions concerning two or more defects. It was 
found that two equally slow sites deep in the bulk "interact" [l^l , in the sense that the overall current is 
significantly suppressed when they are located near each other. Much less has been investigated when the 
two defects are associated with different rates. Both mean-field methods can easily be extended to study 
such issues, especially when the two sites are near each other. At the other extreme, we face a completely 
inhomogeneous TASEP. But this is precisely the scenario more relevant for protein synthesis in vivo. In 
this sense, there is much to be done before we reach the goal of a realistic model for understanding the 
biological process of translation. 
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